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Abstract 

We carry out simulations of a lattice London superconductor in a magnetic 
field B, with a finite magnetic penetration length A. We find that supercon- 
ducting coherence parallel to B persists into the vortex line liquid. We argue 

that the length scale relevant to this effect is A = (/>o/87rT m . 
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In high T c superconductors, thermal fluctuations are believed to melt the ground state 
vortex line lattice at temperatures well below the mean field H c2 line [0,0. The resulting 
vortex line liquid has received intense theoretical and experimental study. In particular, 
recent "flux transformer" experiments |3|.|4]] on YBCO show that superconducting coherence 
parallel to the applied magnetic field B exists over very long length scales well into the 
vortex liquid. 

To investigate the fluctuating vortex line system, numerical simulations have been carried 
out. These have either involved a simplified model of vortex line interactions ||, or have 
used the approximation |§-|8| that the bare magnetic penetration length A — > oo, so that the 
magnetic field B inside the superconductor is uniform. Such A — > oo simulations [§] show 
a sharp transition within the vortex line liquid, corresponding to the onset of coherence 
parallel to B. While this A — > oo model is suggested by the large values of k = A/£ in the 
high T c materials, it may fail close to T c , where the correlation length may exceed the finite 
bare A; thus the finite A model may display different critical behavior from the A — > oo limit. 
To investigate this possibility, we present here new simulations of a system of fluctuating 
vortex lines, in which we include the effect of magnetic screening on the vortex interactions 
||, due to a finite A. We investigate the presence of superconducting coherence within the 
vortex line liquid, and discuss the length scale relevant for this effect. 



Our model is a discretized lattice superconductor in the London limit ||10|| . For simplicity, 
we consider isotropic couplings. Following Carneiro et al. @, a duality transformation maps 
this model onto one of interacting vortex lines, with Hamiltonian 

H = 27r 2 J A 2 5>(r l ) ■ n(r i )G(r i - r,). (1) 

y 

Here n a {vi) (a = x, y, z) is the integer vorticity through plaquette a at site r$ of a cubic 
mesh of points, and G(r) is the lattice London interaction, with Fourier transform 

G g = l/(1 + A 2 Q 2 ), Q 2 = ]T [2-2cos?,J. (2) 

fi=x,y,z 

The coupling is Jo = 0o^o/(16vr 3 A 2 ), with O the flux quantum, and £o the spacing of the 
discrete mesh, which we identify with the bare coherence length that determines the size of 



a vortex core. Henceforth we measure lengths in units of £o, and temperature in units of Jo- 
Periodic boundary conditions in all directions are taken. 

To test for superconducting coherence, we consider the helicity modulus T /J (g Jy ) (where 
fj, ^ v), defined as the linear response coefficient giving the supercurrent j induced by 
a transverse perturbation in the vector potential of the externally applied magnetic field, 

JM = -rM5A?{q v ). (3) 



We have earlier derived JTT| an expression for T^(q^) in a continuum London supercon- 



ductor. The generalization to the lattice superconductor is 

_ [J A 2 ]Q 2 f An 2 J \ 2 {n a {q v )n a {-q v )) \ 
^ 1 + A 2 Q 2 \ VT 1 + A 2 Q 2 J ' [ ) 

where /i, u, a are a cyclic permutation of x, y, z, and Q 2 = 2 — 2 cosq u as g M , q a = 0. 

In the absence of vortices, T M (g y ) = [J A 2 ]Q 2 /(1 + A 2 Q 2 ) describes the total screening of 
the Meissner state. With vortices, an expansion in powers of Q 2 , 

{n a {q v )n a {q v )) = n + n x Q 2 + n 2 Q 4 + (5) 

leads to an expression for T^(q u ) at small q u , in terms of a renormalized coupling [JA 2 ] m r 
and penetration length A^, 

T (a ) = ^ JX2 ^ R (Q) 

with 

[J^ 2 U , 47T 2 J A 2 
Kr i , Ar[2 J o ( n i - n oA 2 ) . . 

+ W — — ' () 

To see the meaning of 7^ and X^r, note that the current j^qv) produced by 
5A^ t (q u ), induces a magnetic vector potential due to Ampere's Law, which in our units 
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is [J X 2 }Q 2 5A™ d (q u ) = j^(q u )- The change in the total magnetic field due to the external 
perturbation is therefore given by 5A t ° t (q u ) = 5A™ d (q v ) + 5A^ !t (q u ), or 

Q 2 



( 9 ) 



Q 2 + Kr\ 

Thus a fraction 1 — 7^ of 5A e ^ t {q u ) penetrates into the material; the remaining fraction 7 M is 
screened out, on the length scale X^r. Equivalently, X^r is the length on which fluctuations 
in magnetic field decay to equilibrium. For a perfect Meissner effect, 7^ = 1, and X^r agrees 



with the usual definition of the London penetration length fl2|l . We therefore interpret 
1/X 2 R ~ p SM as the density of superconducting electron pairs in direction /t, even in the 
more general case of a partial Meissner effect in the mixed state. Although T M (g y ) will 
have the same form Eq.(^) in both the superconducting and the normal metal state (with 
ordinary fluctuation diamagnetism) , a transition will be signaled some singularity in 7^ and 
Xfj,R- We focus now on 7^. 

Consider a uniform applied H = Hz. T x (q y ) and T y (q z ) then describe the response to 
external fields 5H z (q y ) and 5H x (q z ), which represent compression and tilt perturbations of 



H, respectively. Correspondingly one finds [111 



\- lx = dB z /dH z and 1 - ly = dB x /dH x , (10) 

where these susceptibilities are evaluated at the applied field H. Since the high T c materials 
display strong fluctuation diamagnetism even in the normal state, it is unclear whether or 
not j XjV will display a pronounced feature at the superconducting transition. 

For behavior along z however, parallel to H, the criteria for a superconducting transi- 
tion is more clearly defined. T z (q x ), describes the response to an external field SH y (q x ), 
representing a combined shear and tilt perturbation of the uniform applied Hi. In Ref. JTTJ] 
we showed that for the vortex line lattice one has 7 2 = 1 and a perfect Meissner screening. 
For a normal vortex line liquid however, 7 2 = 1 — dB x /dH x 1. Thus the transition to 
the normal state is signaled by a discontinuous jump in 7 2 from unity. Expressed in terms 
of the expansion coefficients of Eq.(^), we have superconducting coherence along the field 
provided n$ = 0; for the normal state, hq > 0. 



For our Monte Carlo simulation, we start with a fixed density / = -B/00 of straight 
vortex lines parallel to z, giving the ground state configuration for an internal field Bz. 
Following Carneiro et al. Q/j, we update the system, heating from the ground state, by adding 
elementary closed vortex rings (a square ring of unit area) with random orientation and 
position. These excitations are accepted or rejected according to the standard Metropolis 
algorithm. This provides a complete sampling of phase space for the vortex variables n(rj), 
subject to the constraints that vorticity is locally conserved, and the average internal field, 
B = (<f>o/V) Ei n(r.j) = f(f) z, is constant. 

Our simulations are for the case / = 1/15, whose ground state on a cubic mesh is a 
close approximation to a perfect triangular lattice. We choose A = 5, comparable to the 
vortex line spacing, a v = l/v7 — 3.87. We study system sizes L± = 30 in the xy plane, 
and L z = 15 and 30 parallel to H. Each data point is typically the result of 5, 000 sweeps 
to equilibrate, followed by 8 — 16, 000 sweeps to compute averages, where each sweep refers 
to L\L Z attempts to add an elementary vortex ring. 

In Fig.[T] we show a sample of our data, plotting {n y (q x )n y (—q x )) vs. q x , for various T, 
and L z = 30. Fitting to Eq.(|5|) through 0(Q 4 ) yields the solid curves, and determines the 
parameters no and n\. Eq.(^) then gives the couplings 7^, which we plot vs. T in Fig.|2|. 
We see that j Xjy decrease towards zero at T m ~ 1.2, while 7 2 decreases at T c ~ 2.0. We also 
show our results for L z = 15. 

In Fig.|3] we show intensity plots of vortex correlations within the same plane perpendic- 



where q_i_ = (q x ,q y ), for various T, and L z = 30. Below T m we see sharp Bragg peaks of 
a vortex line lattice. Above T m we see behavior characteristic of a liquid. T m is thus the 
melting transition. Since the discrete mesh of our simulation acts like a periodic pinning 
potential for vortices, T m also coincides with a depinning of the vortex lines. We believe 
that the drop in 7^ at T m is more a result of this depinning, rather than a direct result of 



ular to B 




(11) 



5 



melting. We expect from Eq.fllOD that j X}y = 1 — dB x /dH x is finite above T m , but this value 
is too small for us to determine accurately. 

With respect to coherence along H, we expect a discontinuous jump in 7 2 from unity to 
1 — dB x /dH x pd at T c . The finite width of the decrease observed in our data is a finite size 
effect; we see that the transition sharpens as L z increases. We therefore estimate T c ~ 2.0, 
well into the vortex line liquid. This is the main result of our simulations. 

Recent flux transformer experiments on YBCO HQ show that there is a temperature 
u T t h" below which vortex line correlations parallel to H become comparable to the thickness 
of the sample. "T^" is clearly above the "Tj rr " where resistivity transverse to H vanishes. 
Resistivity parallel to H however appears 0] to vanish at "T^". A similar conclusion con- 
cerning vortex line correlations may be inferred from the measurements of Ref. [ 13| , where 



the onset of pinning by twin grain boundaries is shown to occur distinctly above a sharp 
first order melting transition. If we identify "T^" with our T c , and "Tj rr " with our T m , our 
results are in complete accord with these experimental findings. 

We have also tried to compute the lengths A^ R , using Eq.(|8|) and our fitted n and n\. 
However the factor [n\ — UqX 2 )/^^ that appears in Eq.(^) is, in the region of the transition, 
the quotient of two small numbers each with large relative error. We were therefore unable 
to obtain meaningful results for X^r. 

Although we have simulated with L z ^> A, one can still question whether our results 
represent the true thermodynamic limit. In particular, in Refs. |3|,§j], it was found that "T^" 
decreased towards "T irr " as sample thickness L z increased. This is also consistent with our 
results for 7 2 in Fig.|2|, where there is some suggestion that, in addition to a sharpening of 
the transition as L z increases, T c also decreases. It is therefore important to note that there 
is another length in the problem 

A(T) = 0o/(8ttT) = (2tv 2 J /T)k\. (12) 

For our simulation, we have 30 = L z A(T m ) ~ 410. 



It has been argued |L4| that A determines the length on which phase correlations C(r) 



^ e »[ ( r )- (o)]^ decay in the vortex line lattice; however these same calculations show that in 
the vortex liquid, the decay length of C(r) is comparable to the spacing between vortex lines 
a v <C L z . Thus this analysis of C(r) does not indicate why A should be an important length 
above T m , where we continue to see superconducting coherence. 

Another possibility is suggested by Nelson's analogy between vortex lines and the 
imaginary time world lines of 2D bosons. Nelson argued that there should be a Kosterlitz- 
Thouless (KT) superffuid transition of the analog bosons. For L z sufficiently large, this KT 
transition would occur at a T c < T m , and so would be pre-empted by the formation of the 
vortex line lattice. But for L z small enough, T c > T m , and one has a new state intermediate 
between the vortex line lattice and the normal vortex line liquid [|15|]. We now restate our 
earlier calculation JTTJ of this T c in order to show that the length which distinguishes between 
these two possibilities is A(T m ). 

As shown by Pollock and Ceperley the 2D boson superfluid density can be expressed 
in terms of the "winding number" W which is the net spatial distance traveled by the 
ensemble of bosons as they travel down the time axis of their world lines. One has p h ° son = 
r mTboson{W 2 ) /2h 2 . According to the KT theory, the 2D superfluid transition occurs at a 
universal value of p b s oson , which translates into the condition (W 2 ) = 4/tt. In terms of vortex 
lines, W just measures the net vorticity transverse to H |L1| . We thus have [T7] 

{W 2 y ) = \im-L(n y (q x )n y (-q x )) (13) 

This is precisely the same correlation as enters T^, and comparing with Eq.fl7|) we can write 
7, = 1 - (A(T)/L Z )(W 2 ). Thus p h ° son oc (W 2 ) = implies 7 2 = 1; the normal boson fluid 
corresponds to a vortex line liquid with coherence along H. p b ° son > implies 7 2 < 1; the 
boson superfluid corresponds to the normal vortex line liquid [[18|j . 

We have shown [JXXJ] that for a normal vortex line liquid, (n y (q x )n y (—q x )) = 
f 2 L\L z T I 'c i4: (q x ), where c 44 (0) = (B 2 /4:7i)dH x /dB x is the tilt modulus, and / = B/(f) . 
Thus we conclude that (W 2 ) = (W 2 + W 2 ) = 2(W 2 } = (8nL z T/(j) 2 )dB x /dH x . This yields 

rp = <Pl dH x n = 4 A(T m ) dH x 

° 2tt 2 L 7 dBr ' Tm ix L 7 dB r 



For large B ^> H c i, dH x /dB x ~ 1. Thus for L z < (4/7r)A(T m ), one has T c > T m and 
hence a vortex line liquid with superconducting coherence along H, intermediate between 
the vortex line lattice, and the normal vortex line liquid. Only for L z > (4/ir)A(T m ) will this 
intermediate state disappear |I9|| . For YBCO, with T m ~ 90° K, one has A(T m ) ~ 1400/im, 
much thicker than the samples (~ 50yum) in Refs. P,P||i~5|. 

We note that for B ^> H&, A(T m ) is a factor B/H c \ larger than the "entanglement" 
length originally proposed by Nelson fl[] as the criterion for the 2D boson superfluid tran- 
sition. This is because the notion of superfluidity, as measured by W, does not precisely 
correspond to the geometric notion of line entanglement. If just as many lines wander to 
the right as wander to the left, one has W = 0, although the lines may still be quite twisted 
and geometrically entangled. 

Finally, we note that for our model, Eqs. (|l^ , |i4"D would predict T c = (8irJ K\/L z )T m ~ 25, 
much higher than the observed T c ~ 2.0. We believe that this results from a breakdown 
of the boson analogy near our observed T c , due to the proliferation of thermally excited 
closed vortex rings, and intersections between vortex lines, such as cause the transition in 
a B = model. This is clearly the case for the A — > oo model where A — > oo at fixed 
Jo, also means A — > oo. Eq.(|TJ]) would then imply T c — > oo. The finite T c found in A — ► oo 
simulations || therefore indicates such a breakdown of the boson analogy. Our present 
results concerning vortex ring distributions and line intersections are similar to what was 
found for A — > oo 0. 

To conclude, we find from simulations with A<I 2 < A(T m ) (as is the case for recent 
experiments) that superconducting coherence parallel to H persists into the vortex line liquid 
state. We argue that this effect should vanish once A(T m ) < L z . 

We would like to thank A. P. Young for very helpful discussions. This work has been 
supported by DOE grant DE-FG02-89ER14017. 
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FIGURES 

FIG. 1. (n x (q y )n x (—q y )} / ' L\ vs. q y = 2irm/L±, (m integer) for various T and L± = L z = 30. 
Sample error bars are shown. Solid lines are a fit to Eq.(|5|). 

FIG. 2. 7^ vs. T for L± = 30 and L z = 15, 30. Sample error bars are shown. ^y x y decreases at 
T m ~ 1.2; j z decreases at T c ~ 2.0. 

FIG. 3. Intensity plot of <S l (qj_) for several T and L± = L z = 30. (a) T = 0.50 < T m shows a 
lattice of Bragg peaks; (b) T = 1.25 w T m ; (c) T m < T = 1.60 < T c in the vortex line liquid; (d) 
T = 2.0 > r c . 
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